On quantum and classical space-bounded processes with 
algebraic transition amplitudes 



John Watrous* 
February 1, 2008 

Abstract 

We define a class of stocliastic processes based on evolutions and measurements of quantum 
systems, and consider the complexity of predicting their long-term behavior. It is shown that 
a very general class of decision problems regarding these stochastic processes can be efficiently 
solved classically in the space-bounded case. The following corollaries are implied by our main 
result for any space-constructible space bound s satisfying s{n) — Q{\ogn). 

• Any space 0{s) uniform family of quantum circuits acting on s qubits and consisting of unitary 
gates and measurement gates defined in a typical way by matrices of algebraic numbers can 
be simulated by an unbounded error space 0(s) ordinary (i.e., fair-coin flipping) probabilistic 
Turing machine, and hence by space 0{s) uniform classical (deterministic) circuits of depth 
O(s^) and size 2*^^*^. The quantum circuits are not required to operate with bounded error 
and may have depth exponential in s. 

• Any (unbounded error) quantum Turing machine running in space s, having arbitrary algebraic 
transition amplitudes, allowing unrestricted measurements during its computation, and having 
no restrictions on running time can be simulated by an unbounded error space 0{s) ordinary 
probabilistic Turing machine, and hence deterministically in space O(s^). 

We also obtain the following classical result: 

• Any unbounded error probabilistic Turing machine running in space s that allows algebraic 
probabilities and algebraic cut-point can be simulated by a space 0{s) ordinary probabilistic 
Turing machine with cut-point 1/2. 

Our technique for handling algebraic numbers in the above simulations may be of independent 
interest. It is shown that any real algebraic number can be accurately approximated by a ratio 
of GapL functions. 



1 Introduction 

In this paper we investigate the complexity of predicting the long-term behavior of stochastic pro- 
cesses induced by evolutions and measurements of discrete quantum mechanical systems. The 
processes considered, which we call selective quantum processes, describe the classical outputs ob- 
tained when operations called selective quantum operations are iterated on finite state quantum 
systems. Quantum Turing machine and quantum circuit computations may be viewed as specific 
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examples of such processes. Selective quantum operations are quite general and include unitary 
evolutions and positive operator valued measures (POVMs). The main result proved in this paper 
regards the space required for classical machines to solve decision problems based on selective quan- 
tum processes; we postpone the formal statement of this result until after necessary background 
material has been discussed. 

Although quantum computation offers the potential for exponential speed-up over classical com- 
putation for certain problems, such as integer factoring and discrete logarithms ||2^, the analogous 
situation does not hold with regard to the amount of space required by quantum machines vs. clas- 
sical machines. It was proved in [24| that quantum Turing machines satisfying certain restrictions 
(discussed below) and running in a given space bound s (for s space-constructible and satisfying 
s{n) = O(logn)) can be simulated by space 0{s) probabilistic Turing machines in the unbounded 
error setting. The assumptions made on the quantum machines were that only rational transition 
amplitudes could be used, and only a very limited class of observations during the computations 
were permitted. By applying our main result to quantum Turing machine computations, we extend 
the above result to machines allowing algebraic transition amplitudes rather than just rational ones, 
and further to machines allowing unrestricted measurements during their computations. Our main 
result also has implications for bounded-width quantum circuits: any space 0{s) uniform family of 
quantum circuit acting on s qubits that consist of unitary gates and measurement gates defined by 
matrices of algebraic numbers (as described later in Section can be simulated by an unbounded 
error space 0{s) probabilistic Turing machine. 

A well-known result of Borodin, Cook, and Pippenger states that any unbounded-error 
space s probabilistic Turing machine computation can be simulated deterministically by uniformly 
generated depth O(s^) circuits with size polynomial in 2*, and hence in deterministic space O(s^) 
(i.e., PrSPACE(s) C NC'^{2'') C DSPACE(s2)). Thus, our results imply the existence of efficient 
parallel (classical) algorithms (in the sense that NC represents a class of efficiently parallelizable 
problems) for predicting the long-term behavior of quantum systems of modest size. It is interesting 
to note that the complexity of the algorithm implied by our technique is independent of the running 
time of the quantum process: in (parallel) time proportional to s^, we may predict the behavior of 
a quantum system consisting of s components (e.g., qubits) that runs for an arbitrary number of 
steps. 

The technique we use to prove our main result is similar to one used in |24|, and has previously 
been used to prove results in classical space-bounded computation (for instance in 1^, |^, |^]). Essen- 
tially, the technique is to manipulate matrices that govern the stochastic processes being considered 
in order to predict their long-term behavior (rather than simulating the processes directly), with 
the matrix manipulations being performed in a very space-efficient manner. In the present case, 
we must modify this technique in order to handle algebraic matrices rather than rational ones, and 
we must reformulate the quantum processes we are considering in the framework covered by the 
technique. Section ^ describes this in detail. 



2 Quantum processes 

In this section we briefly review certain facts from quantum computation and state the definition 
of selective quantum processes that will be used throughout this paper. For a more thorough 
treatment of quantum computing, we refer the reader to the surveys of Berthiaume |^] and Kitaev 
|16|, and to the references therein. Our definition of selective quantum operations is implicit in 
||7|l^. A number of the claims made in this section have straightforward proofs using matrix 
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analysis (see, e.g., Horn and Johnson |1S]), which we omit. 

We restrict our attention to quantum systems having finite classical state sets; for a given 
system, we generally denote the classical state set by S. For example, in the case of quantum 
circuits the set S may be the set of all 0-1 assignments to the wires at some particular level in the 
circuit, while in the case of quantum Turing machines S may be the set of all configurations of the 
machine subject to some given space bound. Given a quantum system with fixed classical state 
set S, a pure state (or superposition) of the system is unit vector in the Hilbert space £2(3). We 
use the Dirac notation to represent elements of ^2('S'); for each s £ S, \s) represents the unit vector 
corresponding to the map that takes s to 1 and each s' ^ s to 0. Elements of ^2(5') are generally 
denoted {ip), \(j)), etc., and may be specified by linear combinations of elements in the orthonormal 
basis {\s) : s G S}. Corresponding to each \^) is a linear functional that maps each vector \(f)) 
to the inner product (conjugate- linear in the first argument). 

A mixed state is a state that may be described by a distribution on (not necessarily orthogonal) 
pure states. Intuitively, a mixed state represents the quantum state of a system given that we have 
limited knowledge of this state. A collection {{pk, \fpk))} such that < pk, YlkPk = 1) ^'^'^ each 1^/;^) 
is a pure state is called a mixture: for each /c, the system is in superposition \tjjk) with probability p^. 
It is the case that different mixtures may yield identical states, in the sense that no measurement can 
distinguish the mixtures even in a statistical sense. For a given mixture {(pfc, \ ipk))}^ we associate 
an l^l X l^l density matrix p having operator representation p = X]fePfc|^fc)(^fc|- Two mixtures yield 
different density matrices if and only if there exists a measurement that can statistically distinguish 
the two mixtures, and so we interpret a given density matrix p as being a canonical representation 
of a given mixed state. Necessary and sufficient conditions for a given l^l x IS"! matrix p to be a 
density matrix (i.e., represent some mixed state) are (i) p must be positive semidefinite, and (ii) p 
must have unit trace. 

A selective quantum operation is a probabilistic mapping that takes as input a density matrix 
p and outputs a collection of pairs (i,p^*-*), each with some probability p^; each is a density 
matrix and i is a classical output that we take to be an integer for simplicity. The output i may be 
the result of some measurement, although this is not the most general situation (for example, the 
system may be measured and part of outcome may be discarded). A selective quantum operation 
£ must be described by a collection {Ai^j | < i < m, 1 < j < /} of IS*] x \S\ matrices satisfying 
the constraint X^^^Lq Si=i ^1 ~ Given such a collection of matrices, we define a function 
Pi : C"^" ^ [0, 1] and a partial function Ei : C"^" ^ C"^" as follows: 

I 



(In case Pi{p) = 0, Ei{p) is undefined.) Now, on input p, the output of £ is defined to be {i,Ei{p)) 
with probability Pi{p) for each i. It may be verified that for any density matrix p, we have < Pi{p) 
and YliPiip) = 1, and furthermore that each Ei{p) is a density matrix (and so the above definition 
is sensible). We also define functions Fq, . . . , Fm as Fi{p) = X]j=i ^i,jP^\ j- '^i^^ simplify matters 
when calculating unconditional probabilities to consider these functions. 

Finally, a selective quantum process is a stochastic process {Rt \ t G N}, where each Rt is a 
random variable whose value corresponds to the classical output of a selective quantum operation. 
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A selective quantum operation £ = {Aij \ < i < m, 1 < j < 1} and an initial density matrix 
Pinit induce a selective quantum process {Rt\t G N} as follows: for n > 1 and ri,... ,r„ G 
{0, . . . , m}, the probability that Ri, . . . , i?„ take values ri, . . . , r„ is the probability that, if the 
selective quantum operation £ is iterated n times given initial state pinit, the resulting classical 
outputs will be ri, . . . ,r„. Thus, the probability that i?„ takes a particular value r„ depends on 
the values taken by . . . , Rn-i in the following way: 

Pr [Rn = rn\Rl = n, . . . , Rn-l = r„_i] = p^,, {Er^^i 0---0 Er-, {pinit)) ■ 

It may be proved by induction that for a selective quantum process {i?^ 1 1 G N} induced by £ and 
Pinit, we have Pr[i?i = n, . . . , i?„ = r„] = tr(Fr„ o • • • o {pinit}) for ah n > 1. 




3 GapL functions and PL 

Next, we recall some definitions and facts from counting complexity and space-bounded complexity. 
Counting complexity is a powerful technique that has its origins in the work of Valiant 
has had a number of applications in complexity theory (including in quantum computing 
For further information on counting complexity, see the survey of Fortnow and the references 
therein. Counting complexity was applied to space-bounded computation in to which the 

reader is referred to for proofs of the theorems stated in this section. For more general background 
information on space-bounded computation, see Saks [f^] . 

Consider a nondeterministic Turing machine M running in logspace. On each input x there 
are some number of computation paths that lead to an accepting configuration and some number 
of paths that lead to a rejecting configuration; we denote these numbers by #M(x) and 7^M(x), 
respectively. 

Definition 3.1 A function f : T,* ^ I, is a GapL function (f G GapL^ if there exists a logspace 
nondeterministic Turing machine Mj such that f{x) = ^Mf{x) — ^Mf{x) for every input x. 

GapL functions characterize the class PL as follows: 

Theorem 3.1 Let yl C S*. Then A G PL if and only if there exists f G GapL such that x & A 4^ 
/(x) > for every x G H*. 

Any integer function computable in logspace is necessarily a GapL function: FL C GapL. 
GapL functions are quite useful for proving results in space-bounded complexity due to the fact 
that various closure properties of GapL functions hold: 

Theorem 3.2 Let f G GapL, let g G FL, and let p be an integer polynomial. Define functions hi, 
h2, and h^ as follows: hi{x) = f{g{x)), h2{x) = X]f=o'^ /(^' 0; h^ix) = nf=o'^ /(^' ^^en 
hi, /i2, and h^ are in GapL. 

Here it is assumed that nonnegative integers are identified with their encodings as binary strings 
in the usual way. 

There is a close relationship between GapL functions and the determinant function, which is 
crucial for the proof of our main theorem: 

Theorem 3.3 Let f G GapL and let p he any integer polynomial satisfying p[n) > 1 for rz > 0. 
For each x G S* let A{x) he a p{\x\) x p(|x|) matrix defined as A{x)[i^j] = f{x,i,j), and define 
g{x) = det(^(x)). Then g G GapL. 

See Q for a proof of this theorem. 
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4 Main result and applications 



Now we are prepared to state the main theorem. 

Theorem 4.1 Let p and q be integer polynomials satisfying p{n),q{n) > 1 for n > 0, let 
{ai,... ,afc,/3} be a finite collection of algebraic numbers with (3 real and in the range [0,1], and 
let ai, bi, ci, and di (for 1 < I < k) be GapL functions such that each bi and di is nonzero on 
all inputs. For each x £ S*, let be a selective quantum operation described by p{\x\) x p(|x|) 
matrices {A^^ij | < i < 1 < j < gdxl)} defined as follows: 

k , . . ./ .f\ 

and let px be a p{\x\) x density matrix specified as follows: 

r-/ -n ciix,i',j') 

Let {Rx,t I i € N} &e the selective quantum process induced by Ex and px for each x. Then the 
language 

{x I Pr[3 1 : Rx,i = ■■■ = Rx,t-i = 0, Rx,t = 1] > /?} 

is in PL. 

This theorem is sufficiently general to imply that a wide variety of "logspace quantum computa- 
tions" can be simulated in PL, including logspace quantum Turing machine and logarithmic-width 
quantum circuit computations that have algebraic transition amplitudes and allow measurements 
during their computations. We now discuss these special cases in more detail. We conclude this 
section with a brief discussion on how these results can be extended to a more general class of 
space-bounds. 

We will begin with quantum circuits, since it will be straightforward to translate results in this 
setting to the quantum Turing machine model. We use the quantum circuit formalization given in 
lH] , which allows a very general class of quantum gates including unitary gates and "measurement 
gates". In order to fix our notation, we first briefly review this formalism. 

A k qubit quantum gate is a linear mapping on the space of 2^ x 2^ complex matrices induced by 
a collection of 2^ x 2^ matrices {Ai, . . . , A^} in a similar manner to selective quantum operations 
(but with no classical output). Specifically, we require Ylj^]^j — ™d have that this 

collection induces the mapping p i— > Ylj ^jP^] density matrices. A k qubit gate acts on an 
ordered /c-tuple of qubits in the natural way: if we associate the initial mixed state of the qubits 
with a 2^= X 2^^ density matrix p, the effect of the gate is given by the above mapping. An n qubit 
quantum circuit is a sequence of quantum gates applied to ordered subsets of a collection of n 
qubits. The action of a A; qubit gate on a given ordered subset of n qubits, for n > A;, is given by 
taking the Kronecker product of each matrix Aj with the 2*^"^ x 2""'^ identity matrix, permuting 
rows and columns of the resulting matrices according to which qubits are acted on by the gate, and 
applying the resulting operation to the state of the n qubits. 

We now define what we mean by a logspace-uniform family of quantum circuits acting on a 
logarithmic number of qubits. Let Q = {Gi, . . . , Gk} be a fixed, finite collection of quantum gates. 
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each acting on a constant number of qubits and specified by a collection of matrices as above. We 
require that each such matrix have entries that are algebraic numbers. Let s be a space-constructible 
function satisfying s(n) = 0(logn). Since s is sub-linear, we assume the particular input x is not 
given as input to the quantum circuit, but rather is input to the deterministic procedure that 
generates the circuits. The input to the circuit is assumed to be s{\x\) qubits each in the |0) state. 
Let / be a mapping that takes an input x and an integer in the range {1,... ,r(|x|)} as input, 
for some polynomial r, and outputs an index of an element in Q along with an ordered subset of 
{1,... ,s(|x|)}, specifying a gate and the qubits upon which that gate is to act. The quantum 
circuit generated by / is the sequence of gates given by f{x, 1), . . . , f{x,r{\x\)), applied in order. 
One of the qubits is specified as the output of the circuit, and is assumed to be observed in the 
0-1 basis after the circuit has been applied, yielding acceptance or rejection. (In case we wish to 
consider the output of the circuit to be a function, multiple qubits may be specified as output 
qubits.) If the function / can be computed in logspace by a deterministic Turing machine, then we 
say the resulting family of circuits is logspace-uniform. 

Now, we claim that for any logspace-uniform family of quantum circuits acting on a logarithmic 
number of qubits, the language consisting of those inputs x accepted with probability exceeding a 



given algebraic cut-point /? is in PL, following from Theorem 4.1. This may be shown by defining 



a selective quantum operation acting on s{\x\) + [log2r(|x|)] qubits: the first s{\x\) qubits 
represent the qubits in the quantum circuit, and the remaining qubits index the gates in the 
circuit. The selective quantum operation £x effectively measures the qubits indexing the particular 
gate to be applied, applies that gate appropriately to the first s qubits, increments the index 
qubits modulo and outputs one of the classical results (the computation has not yet 

completed), 1 (the circuit accepts), or 2 (the circuit rejects). Given that the mapping / is logspace 
computable, it is possible to define GapL functions a/ and bi (in fact, each 5/ may be constant, 
taking value 1, and each a/ may be an FL function) such that £x is given by p(|x|) x pd^l) matrices 
{^x,i,j I < i < (7(|x|), 1 < j < ^d^l)} for polynomials p and q as in the statement of Theorem [4.1| . 
(Since the output of £x is always 0,1, or 2, we will have A^^tj = for i > 2.) The probability 
that the circuit accepts is precisely Pr[3 1 : R^^i = • • • = Rx,t-i = 0, Rx,t = l], for {Rx,i-,Rx,2-, • • • } 
the process induced by £x (along with the initial density matrix px describing the initial zero state 
of the qubits and initial state of the index qubits — easily seen to be computable in the sense of 



Theorem |4.1| ) . A more formal presentation of this construction will appear in the final version of 
this paper. 

Next we discuss quantum Turing machine computations. We are not aware of any systematic 
treatment of quantum Turing machines that may perform unitary operations and measurements 
during their computations, nor will we attempt to provide such a treatment here. However, we 
claim that any reasonable notion of a quantum Turing machine M running in logspace that allows 
measurements during its computation may be formulated in terms of a selective quantum process 
in such a way that the following holds: the language consisting of all strings accepted by M 
with probability exceeding some algebraic cut-point (3 reduces to the language in the statement of 



Theorem |4.1| , and hence is contained in PL. 

For instance, consider a Turing machine consisting of two parts: a classical part and a quantum 
part. The input tape may be considered to belong to the classical part, along with a classical work 
tape and a classical portion of the internal state, while the quantum part consists of a quantum 
work tape and a quantum portion of the internal state. For each local description of the classical 
part of the machine, we may have a quantum transition function that specifies the evolution of the 
quantum part of the machine in the usual manner (e.g., as described by Bernstein and Vazirani pf). 
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Such "quantum steps" may be alternated with "classical steps" , in which the classical part of the 
machine evolves classically, perhaps involving measurements of the quantum portion of the internal 
state. Under the assumption that the specifications of the quantum transition functions and the 
measurements of the quantum portion of the internal state are described by finite collections of 
algebraic numbers, the problem of determining if such a machine accepts with probability exceeding 
some cut-point (j reduces to the problem in Theorem 4.1. This may be argued by referring to the 
previous discussion on quantum circuits. Given such a quantum Turing machine, we may represent 
the work tapes, internal state, and input tape head position of this machine by the qubits in a 
logspace uniform quantum circuit, with the read-only input of the Turing machine corresponding to 
the input to the logspace function generating the circuit. Similar to the quantum circuit simulation 
of quantum Turing machines due to Yao [25|, we may define a quantum circuit that simulates 
one step in the Turing machine computation. As above, we may then define a selective quantum 
operation £x that simulates the action of this circuit, and produces classical output 0, 1, or 2 as 
above, with output 1 or 2 corresponding to the situation that the Turing machine has entered an 



accepting or rejecting state, respectively. Here we take advantage of the fact that Theorem 4.1 



places no restriction on the running time of the quantum process, since the Turing machine being 
simulated need not necessarily halt absolutely or even with probability 1. Again we postpone the 
formal presentation of this construction to the final version of this paper. 

We also note that probabilistic Turing machines having algebraic probability transitions and 
algebraic cut-point are a restricted case of the quantum Turing machines we have considered. 
Theorem 4.1 thus implies that even in the unbounded error setting, logspace probabilistic Turing 
machines having algebraic transitions and cut-points are equivalent in power to ordinary (i.e., 
fair-coin flipping) logspace probabilistic Turing machines with cut-point 1/2. 

Finally, we mention that the above results may be extended to more general space bounds by 
standard padding arguments, under the assumption that the space bound s is space-constructible 
and satisfies s{n) = r2(logn). Assume that we have a space 0{s) uniform quantum circuit (de- 
fined analogously to logspace uniform quantum circuits) acting on s qubits, and let A be the 
language defined by the resulting circuits given some algebraic cut-point (3. Define a new language 
A = < xQl X ^ A> . It Is straightforward to show that this language has logspace uniform 

quantum circuits, and hence is in PL. This follows from the fact that the suffix i^*""''^ easily 
be recognized and ignored in logspace (following from the fact that s is space constructible, so we 
may simulate the machine that marks of s(|a;|) tape squares and reject if this simulation requires 
more than logspace), after which the original computation is performed on the prefix x in space s, 
which is logarithmic in the length of x01^° \ Given that A is in PL, it is also straightforward 
to show that A is in PrSPACE(s) by similar arguments. For a more thorough discussion of such 
techniques, see, e.g.. Section 6.4 of pF 



5 Proof of the main theorem 



In this section we present a proof of Theorem ^T 

proof. We then present various facts needed for the formal proof in Sections [5.1| and |5^ 



We begin by outlining the main ideas of the 

and 



assemble the parts in Section |5^ 



The method used to prove the main theorem is similar to one used in a number of other papers 
on space-bounded computation, particularly in Q and, in the quantum setting, [24|; in short, the 
long-term behavior of a given selective quantum process is determined by performing various matrix 
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operations on a matrix that determines the behavior of the process. 

Consider first the simpler case of stochastic processes described by Markov chains. Assume 
the Markov chain has state set {1, . . . ,N}, state 1 is the initial state, and states — 1 and 
are absorbing states. Assume the chain is described by a transition matrix A, and let i? be a 
modification of A where columns A^ — 1 and A^ are set to zero. The probability that the process 
eventually enters state N is given by -B* [A^, 1] . Under the assumption that B has eigenvalues 
strictly less than 1 in absolute value (i.e., the process eventually enters state A^ — 1 or A^ with 
probability 1), the probability that the process eventually enters state A" is given by 

(We write Xij to denote the matrix obtained by removing the ith row and jth. column of a given 
matrix X throughout this paper.) Thus, the problem of determining if this probability is strictly 

larger than 1/2 reduces to determining the sign of (— 1)^+^ '^'^det(J-^)^"^ ~ h- 

In the situation that the chain above depends on some input string x, where we assume A^ is 
polynomial in |x| and entries of A are rational numbers computable in logspace, it is possible to 
determine whether the chain associated with x enters state A^ with probability exceeding 1/2 in 
PL as follows: we define a GapL function that takes the same sign as (— 1)^"*"-^ '^°det(J-Bj^'^ ~ I' 



and apply Theorem 3.1. The fact that such a GapL function exists follows from the properties of 



GapL functions given by Theorems 3.2 and |3^ 

Now, in the case of selective quantum processes defined by matrices of algebraic numbers, the 
situation becomes somewhat more complicated, although the main idea is the same. The first 
issue we must face is the arithmetic with algebraic numbers. We do not know how to approximate 
algebraic numbers to a sufficient degree of accuracy for the above technique to work, supposing that 
this approximation is to take place in deterministic logspace. Instead, we approximate algebraic 



numbers by ratios of GapL functions; Section 5.1 below describes what we mean by a sufficiently 



accurate approximation, and proves that this approximation can be achieved. The second issue is 
that it is not immediate that a given selective quantum process is governed by a single matrix as 



in the case of Markov chains. Indeed this is the case, however, as we note in Section 5.2. Here we 
also demonstrate how the basic technique from above can be extended to the resulting matrices, 
which are not necessarily stochastic and may have eigenvalues on the unit circle. 

5.1 GapL approximable numbers 

We now define a class of numbers that can be efficiently approximated by ratios of GapL functions, 
and then show that this class includes the algebraic real numbers. 

Definition 5.1 Let a S M. We say a is GapL approximable if there exist f,g^ GapL such that 
for all n > we have g{l"') ^ and 



/(I") 



a 



5(1") 

Denote the set of GapL approximable numbers by G. 



< 2~". 



Theorem 5.1 Let a be any real algebraic number. Then a € G. 
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The proof of Theorem 5.1 rehes on the following lemma. 

Lemma 5.2 Let uq and ui be bivariate integer polynomials and let ao and ai be integers. Then 
there exists f £ GapL such that 



/(l",c) 
for n > and c £ {0, 1}. 



ar n = 1 



Proof. We will define a logspace NTM Mj such that /(I", c) = #M/(1", c) - #M/(1", c) satisfies 
the recurrence in the statement of the lemma. 

Write Uc{X,Y) = Ylo<i j<d'^c,i,jX^Y^ . To simplify the presentation of Mj, we define and 

Mc,i,j for c G {0, 1}, < i, j < d to be NTMs that take no input and satisfy #Mc — #Mc = Oc and 
^Mc^ij — ^Mc^ij = Uc^ij- The execution of Mf may be described as follows: 



Input: (l",c). 

Set k = [logn] and s = 1. 

Call P(c). 

If s = 1 then accept, otherwise reject. 

Procedure P(c) 

If /c = 0, simulate Mc and set s = —s if Mc rejects. 
Else 

Guess i, j £ {0, . . . , d}. 

Simulate Mc^ij and set s = — s if Mc^ij rejects. 
Set k = k — 1. 
Repeat i times: Call -P(O). 
Repeat j times: Call P(l). 
Set A: = A; + 1. 
End Procedure P. 



The variables k and s are "global", while i, j, and any auxiliary variables needed by Procedure P 
are "local". Since i, j, and all required auxiliary variables are constant in size, Mf will need to 
store only a constant amount of information for each level of the recursion. As the recursion will 
have depth at most logarithmic in n, Mf requires space O(logn) to implement the recursion. Since 
each of the machines Mc and M^^ij require only constant space, it follows that Mf may be taken 
to run in space O(logn). 

Now let us analyze the computation of Mf. Each execution of Procedure P causes the com- 
putation of Mf to branch along several computation paths, each path having the effect of either 
leaving s unchanged or replacing s with —s. Let r~^{k,c) denote the number of computation paths 
induced by calling P{c) for a given value of k that leave s unchanged, let r~{k, c) denote the number 
of computation paths induced by calling P{c) that result in s being replaced by — s, and define 
r(fc,c) = r+{k,c) - r-{k,c). Note that we have #M/(l",c) - #iW7(l",c) = r([logn],c). Since 
[log[n/2]] = [logn] — 1 for any integer n > 2, it remains to prove that r obeys the recurrence 

( Uc{r{k- 1,0), r{k- 1,1)) k>l 
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In case k = 0, Procedure P{c) induces t^Mc computation paths that do not modify s and 
paths that replace s with —s, and thus r(0,c) = Oc- Now suppose k > 1 and assume that the 
number of paths induced by P{b) that do not change s (replace s with —s) when k is replaced by 
A; — 1 is described by r^{k — 1,6) {r~{k — 1,6), respectively) for each 6. For each pair z, j that 
may be guessed, it may be proved (using the binomial theorem) that the number of computation 
paths induced by the remaining portion of P(c) that have the effect of leaving s unchanged minus 
the number of paths that replace s by —s is given by Uc,ijr(k — 1,0)* r(A; — 1, 1)-'. We therefore 
conclude that 



r{k, c) 



for k > 1, which completes the proof. 



Y,Uc,i,jr{k-l,Oyr{k-l,iy 



Uc{r{k-l,0),r{k-l,l)) 



Proof of Theorem 5^. Clearly we have G G, so consider the case a 7^ 0. Let p{x) = 
Pdx'^ + • • • + be an integer polynomial such that p{a) = 0, and assume without loss of generality 
that p'{a) / 0. 

Lemma will allow us to use Newton's Method to approximate a by GapL functions. We 
have that there exist positive constants ^ and K, where and S^K are at most 1/2, such that for 
xo G (a — + and x^+i = — p{xk)/p'{xk) for A; > 0, the inequality \xk+i — a\ < K\xk — 
is satisfied for all A; > 0. Thus we have \xk — a\ < 2" 

Define 



^'^ for every A; > 0. 



and note that 



uo{x,y) = ^{j-l)pjX^y ^ 

j=0 

d 



uoix,y) _ p{x/y) 
ui{x,y) y p'(x/y)' 



Let ao,ai S Z, oi / 0, be such that \a — ao/ai\ < ^. By Lemma 5^ there exists / G GapL such 
that 

/(1",0) ^ no(/(ir"/2l,0),/(irn/2l^i)) ^ /(irn/21^0) p[ ^/^\lZf) ) 

/(I", 1) ni(/(ir"/2i , 0), /(ir^/21 , 1)) /(irn/21 , 1) 

for n > 2, and /(1,0)//(1, 1) = ao/ai- Consequently 



/(1",0) 



/(l^l) 



a 



< 2-2^'°' "1 < 2- 



for every n > 1. We may now define g,h £ GapL that satisfy 



5(1") 



/(1",0) n>l 
oo n = 



10 



and 



u(,n. _ r n>l 
" 1 ao n = 0, 

so that - a| < 2"" for all n > 0. Thus a G G. ■ 

It is interesting to note that the set G is in fact a subfield of the reals, following from a 
straightforward proof relying on the closure properties of GapL functions. It is also straightforward 
to prove that G contains some transcendental numbers (such as tt, for example) so G properly 
contains the algebraic reals. 



5.2 Quantum processes and matrix problems 



Next we prove that selective quantum processes may be described by transition matrices in a 
manner similar to Markov chains. It will simplify matters to note first that the selective quantum 
operations and density matrices underlying a given selective quantum process may be assumed to 
be real. 

Lemma 5.3 Let {Ri, R2, ■ ■ ■} be a selective quantum, process induced by selective quantum opera- 
tion £ = {Aij} and initial state pmit- Define real matrices {^ij} and p'^^^^ as follows: 

4^.[2i'-l,2/-l] = 
4,.[2i'-l,2/] = ^{A,,\i\j']) 
Ar.[2i',2j'-l\ = -9(A,,,.[i',/]) 



and 



i^pi/-l,2j'-l] = i3f^(p,n^^[^^/]) 
[2i'-l,2f] = ^^{pinit[i',j']) 



init 

,j2i',2/] = ^?ftipinit\i',j']) 



^^{Pinit[i',j']) 



Then j} and p^^^^ also induce the selective quantum process {i?i,i?2, . . . }• 

The proof of this lemma is straightforward. 

Recall that for n x n matrices A and B, the Kronecker product yl i? is an x matrix 
satisfying {A (g) B)[{io - l)n + ii, (jo - l)n + ji] = A[io,jo] B[ii,ji] for 1 < io,ii,jo,ji < n. For 
fixed n, let us also define a mapping vcc from n x n matrices to dimensional (column) vectors 
as vec(^)[(i — l)n + j] = A[i,j] for I < i,j < n. It is easy to prove that for n x n matrices A, B, 
and C we have vec(^SC) = (A (g) C^) vec{B) and tr {A'^B) = vec(A)^ vec(S). 

Lemma 5.4 Let {Ai j | < i < m, 1 < j < /} describe a selective quantum operation, let Pmit 

be 

an initial state, and let {Ri,R2, . . . } be the induced selective quantum process. For given f3 define 
an {n^ + 2) x [n^ + 2) matrix M as follows: 



( 



M = 





vec(pj„it) 

-a 







•A 







vec{EjAlA,,y 0) 
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Then all eigenvalues of M are bounded in absolute value by 1. Furthermore, for each nonnegative 
integer t we have M^+^[n^ + 2, 1] = Pr[i?i = 0, . . . ,Rt = 0,Rt+i = 1]- 



The following lemma is used in the proof of Lemma 5.4. The proof is a modification of the proof 
of Lemma 1 in |21|. 



Lemma 5.5 Let {Aij | < i < m, 1 < j < /} satisfy J2i j A-ljAij = I. Then for each i, J2j A 
Aij has eigenvalues bounded by 1 in absolute value. 



Proof. Let f 7^ and A satisfy i^- Aij0Aij)v = Xv, and let B be the matrix such that vec{B) = v. 
AsB ^0, there exists a uni t vect or IV') such that (V^I^IV") / 0. Define C = {ip\B'f\^)B + {'ilj\B\'ip)B^ 
and write \C\ to denote V C'^C. Note that |C| and |C| + C are positive semidefinite, as C is 
hermitian. Let Fi be as defined in Section |2[ so that Fi{B) = XB and Fi{B'^) = XB^ . Thus, we 
have {il^lF^HlCl + C)\^) = + 2\{i;\B\i;)\'^^{X^) for A; > L Since \C\ + C and |C| are 

positive semidefinite, we have < {iplFj'ilCl+Qlip) < triFj'dCl+C)) < tr(|C| +C), and similarly 
< (V'|i^''(|C|)|V') < tr(|C7|). Consequently, 2|(V'|S|V'>p3ft(A'=) is bounded (independent of k). As 
5|^)|2 ^ 0, this implies |A| < 1. ■ 



Proof of Lemma |5.4| . A straightforward computation shows the following: 

i?"+2[n2 + 2,l] = tr(^A,,,FS'{p^n^t)Al 

= ti{FioF^{p,nit)) 

= Pr[i?i = 0,... = 0,i?^+i = 1]. 



Thus it remains to show that all eigenvalues of B are bounded by 1 in absolute value. By Lemma 5.5 



we have that all eigenvalues of Aq^j^Aq^ are bounded by 1 in absolute value. Since any nonzero 
eigenvalue of B must be an eigenvalue of Ylj ^o,i ® ^o,j) the required fact holds. ■ 



Lemmas 5^ and 5^ will allow us to translate the problem in Theorem LI regarding selective 
quantum processes to an equivalent matrix problem. The next theorem proves that this matrix 
problem is solvable in PL. 

Theorem 5.6 Let p be an integer polynomial satisfying p{n) > 2 for n > 0, let = {ai, . . . , a/.} 
be any finite collection of real algebraic numbers, and letri,... ,rk,si,... ,Sk € GapL such that 
each si is nonzero on all inputs. For each x £ T,* let be a p{\x\) x p{\x\) matrix defined as 

for 1 < i,j < p{\x\). Then if has eigenvalues bounded by 1 in absolute value and the series 
X]j>o -^ib(l^l)' 1] converges, we have 



t>0 



Mi[p(|x|),l] >0 V GPL. 
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The proof of this theorem rehes on a few technical facts, which we now state. First, however, let 
us mention some notation: for any univariate polynomial f{X) = fj^'' or bivariate polynomial 
f{X,Y) = Yli j fi-.j^^^^ write ||/|| to denote maxj{|/j|} or maxij{|/jj|}, respectively. 

Lemma 5.7 Let u and v be polynomials of degree at most d > 1 such that \v{0)\ > 6. Then for 
\z\ < e < 2d\\v \ \ ' ''^^ have \v{z)\ > 5/2 and 



n(0) u{z) 



v{0) 



v{z 



< 



Aed \\u\\ \\v\ 



The proof is straightforward. 

Let us now state a theorem due to Mahler that will be used below (see pages 44-46 of [|l7[| for 
a proof). 

Theorem 5.8 (Mahler) Let f and g he integer polynomials of degree df anddg, respectively, and 
let a satisfy f{a) = and g{a) ^ 0. Then \g{oi)\ is greater than or equal to the quantity 



{df+dg-l)\ 



dg \\g\\<^s-^ (|a|'^/-iH h|a| + 1)' 



Lemma 5.9 For any real algebraic number a there exist positive integer constants Ci and C2 such 
that the following holds. Let g and h be bivariate integer polynomials such that ||g||,||/i|| < 2^ and 



deg((7), deg(/i) < A^, for N >2, and such that lim 



^Tl h(a,z) 



exists. Then 



lim 



g{a,z) 



2T1 h{a, z) 

Furthermore, for ^ and a satisfying ^ < 2~^'^^'^ and \a — a\ < ^^Af+i^ have h{a, 1 — ^) 7^ and 



g{a,z) g{a,l-£,) 



lim , , 
zTi h[a, z) h{a, 1 — C) 



< 2 



-CiAf2 



Proof. First, we note that by Lemma |5.7| there exist positive constants C3 and C4, depending only 
on a, such that for any polynomials p and g satisfying deg(p), deg(g) < A^, < 2^^(1 + |q|)^, 

and \q{0)\ > S for 6 > 0, we have \q{i^)\ > S/2 and 



p(0) p{i^) 



< 



\u\ 2^4 A^' 



whenever < 62-^^^ . (Of course these inequalities are not tight — rather they are chosen to 
simplify arithmetic and notation below.) Furthermore, by Theorem ^.8| there exists a positive 
constant C5, again depending only on a, such that for any polynomial p satisfying deg(p) < N 
and \\p\\ < 2^^ we have |p(a)| > 2~'~"^^ whenever p{a) 7^ 0. Without loss of generality assume 
C5 > 2 + 2|a|. 

Now, define u{x) = g{a, 1 — x) and v{x) = h{a, 1 — x). We may write 



N 



u[x) 



^aj(a) 

j=0 



and 



N 
j=0 



13 



for integer polynomials aj and bj, < j < N, satisfying deg(aj), deg(6j) < and ||aj||, ||fej|| < 2^^. 
Let k = mm{j \ hj{a) / 0}. As hk{a) / 0, we have \hk{a)\ > 2-^''^\ Similarly, \ak{a)\ > 2-^''^^ 
in case ak{a) is nonzero. Define uo{x) = u{x)/x^ and vq{x) = v{x)/x^ . As we assume lim^^i 
exists, we must have aj{a) = for j < k, and hence uq and vq are polynomials. Furthermore, we 



have limzti = ^^4^ = itt^ and ^^r^ = "°!^%^ whenever v(x) ^ 0. Consequently 

h(a,z) vo(0) oi^{a) v(x) vo{x) \ / / J 



u{x) uo(x) 



g{a,z) 



zTi n[a, z) 



> 



22^(l + |a| 



y 2-2C5N^ 



(1) 



whenever the limit is nonzero. Also note the following: deg(tto), deg(i;o) < A^, 

\a\)^, and \vo{0)\ = \bk{a)\ > 2-^'^^\ Thus, for ^ < 2-(^3+C5)iv2 foj^Q^g ^j^^^ |^^(^)| > i 2-C5A^2 
and 



^'o(0) t;o(e) 



<^2 



(C4+2C5)Ar2 



(2) 



Now assume ^ < 2 i'^'i+'^s)^'^ [g fixed, and define r(x) = g{a — x, !—£,) and s{x) = h{a — x, 1— '^). 
We have deg(r), deg(s) < A^, ||r||, ||s|| < 22^(1 + |a|)^, and 



1^0)1 



mi > ^e^2-^^^^ 



Thus, for |a - 5| < ^2^+^ < ^ ^^2-(^3+'^5)A^2 conclude that s{a - a) = h{a, 1 - / and 



r(0) r(Q — a) 



By (D and 



s(0) s(a — a) 
we therefore have 



(C4+2C5)Af2 



g{a,z) g{a,l-C) 



lim , , 
zTi h{a, z) h{a, 1 - 



< 



< 4^2 



no(0) (0 



(3) 



+ 



r(0) r(a — a) 



s(0) s{a — a) 



(4) 



Now, define Ci = [2 C5 + 1] and C2 = Ci + [C3 + C4 + 2 C5 + 3] . By (0) we have 



1™ 1, / ^ 
2T1 n(a, z) 



> 2" 



as 



2 CiN^ < I2 2C5Af^_ Furthermore, for ^ and q satisfying ^ < 2 '^^N^ ^^j^j |q, _ q,| < ^ 



2Ar+l 



we 



have ^ < 2-('^3+C5)Af ^ ^nd thus by (|) we have 



gia,z) g{a,l-^) 



lim 

2T1 /ila, z) h{a, 1 — £,) 



7V2 



as required. 



Proof of Theorem |5.6| . First, we outline briefly the main idea of the proof. (Throughout the proof 
we let p denote as |x| is always the point at which p is evaluated). Under the assumption 
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that Mx has eigenvalues bounded in absolute value by 1 and the series X]t>o ^'^x\P^ 1] converges, we 
have 



^ ^ zn det(I 



det((/ - zM^)i,p) 



zMx) 



(5) 



We approximate the algebraic numbers comprising by ratios of GapL functions, and we ap- 
proximate the limit by substituting for z a quantity very close to 1. Relying on the fact that the 



determinants of matrices defined by GapL functions are also in GapL (Theorem 3.3), our approx- 
imation of (^) will be a ratio of GapL functions. Based on these GapL functions, together with a 
bound on the error of the approximation, we define a GapL function F such that F{x) > if and 
only if X^j>o ^x\p^ 1] > 0. By Theorem 3.1, this suffices to prove the theorem. In the remainder of 



the proof, we define the function F and demonstrate that it is indeed the case that F{x) > if and 
only if Xlt>o ^xIp^ 1] > 0- -l^or convenience we assume |x| > 2, since F may be modified on inputs 
of length and 1 without changing the fact that it is a GapL function 

Let a be a real algebraic number such that Q[a] = . . . , otk]] such an a always exists as 

Q[ai, . . . ,afc] is a finite degree (separable) extension of Q (see, e.g., [14|, page 284). Let d be the 
degree of the minimal polynomial of a and fix positive integers m and B and integer polynomials 
,qk so that a/ = deg(g/) < d, and \\qi\\ < B for 1 < I < k. Define bivariate integer 



Qi 



polynomials wi,... ,Wk as wi{yi,y2) = 2/29^(^1/2/2), and note that deg(w/) < d and \\wi\\ < B for 
1 < I < k. Also note that d, m, B, qi, ■ ■ ■ ,qk and wi, ... , depend only on fi, and not on the 
input X. 
Define 



p p k 

h{x) = m JIJJ Jl s/(x,i, j) 
i=i j=i 1=1 



By Theorem 3^, h G GapL. Let Ex{y) be a p x p matrix defined by 



qi{y)ri{x,i,j) 
msi{x,i,j) 



1=1 

For each Ex{y)[i,j] is an integer polynomial in y. We have Ex{a) = h{x)Mx. Next, define 

ux{y,z) = {-1)^+Ph{x) det{{h{x)I -zEx{y))i,p), 
Vx{y,z) = det{h{x) I - zEx{y)). 

Note that there exists a positive integer constant C such that deg{ux) ■, deg{vx) < and 
ll^^xIMI^^xll < 2l^l . Given that M^: has eigenvalues bounded by 1 in absolute value and ^^>Q M*[p, 1] 
converges for each x, we have 



Ux(a, Z) v-^ . r 

hm ^) ' = yM* p,l . 



t>0 



By Lemma 5.£, there exist positive integer constants Ci and C2 such that Vx {a,l — 2 '^al^'P'^ \ ^ 



^M*b, l]-2- 



-Ci|x|2 



> 2 



-Ci\x 



2C 



(6) 
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and 



< 



2-^11x1 



2C 



(7) 



whenever |a — a| < 2 3C2|a;p<^^ 

By Theorem 5J there exist GapL functions / and g such that — a| < 2~" for each 

n > 0. Define = 3C2|xp*^ and write 5 = g^-^i,( \ x \ )-^ ■ The value 5 will be our approximation 

of Q. Also define = C2|xp'^. The value 1 — 2~^(l^l) will be substituted for z in order to 

approximate the limit. 
Next, define 



1=1 



p p k 

i'=i j>=i i>=i 



for 1 < i,j < p, and let denote the p x p matrix defined by ^^^[i, j] = a(x,i,j). Here we let 
7^ denote the value 1 or depending on whether ^ or = 

respectively. By Theorem U, a G GapL. Note that = {g{V^\'=\^))'^Ej^a). Define 

b{x,i,j) = /j(x)(5(l'^(l^-|)))"2^(l-'l)[i = j] - (2^(1-1) - l)a{x,i,j), 

and let be the p x p matrix defined by Bx[i,j] = b{x, Thus, we have 

Bx = /i(x)(5(r(l^l)))'^2'^(l^l)/- (2^'(l^l) - 1)^^. 



By Theorem 3.2, b G GapL. Next, define 

U{x) = (-l)P+i/i(x)(5(l'^(l^l)))^2'^('^')det((5^)i 
V{x) = det{B,). 



pj^ 



By Theorem 3.2 and Theorem 3.3, U,V £ GapL. Finally, define 

F{x) = 2^il^'l'^C/(x)y(x) - 



By Theorem U, F G GapL. 

It remains to show that for every x G S*, F{x) > if and only if Y2t>o ^xlP^ ^- 
above, it may be verified that 

U{x) = (5(l^(l^l)))'^P2P'^(l^l)n^.(5,l -2^(l^'l)) 

and 

V{x) = 2P'^(I^I) 1-^(5, 1 -2'^(l^l)). 

Thus we have V{x) ^ 0. Furthermore F{x) > if and only if 



V{x) > 



which is equivalent to 



(a,l - 2M|a^l)) 

By (|6|) and the inequality in (^) holds if and only if J2t>o ^xIPj M > required. 



(8) 
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5.3 Completion of the proof 



Now Theorem ^?T| follows in straightforward fashion. Let p, q, {ai,... , 0^,(3} , ai, bi, ci, and 
di, 1 < I < k, and {Rx^i, Rx,2j ■ ■ ■} be as in the statement of the theorem. By Lemma 5.3 we 
may assume {ai, . . . , a^} are real algebraic numbers, since otherwise we modify the ai, bi, ci, and 
di functions and take the real and imaginary parts of {oi,... ,0^} accordingly. Let be the 



+ 2) X + 2) matrix described in Lemma |5.4| for each input x. Using Theorem 3.2 it 

is routine to define GapL functions r/ and si, 1 < I < k + 1, such that 



fc+i 



^ si{x,i,j) 



where we write a^+i = /3- By Lemma 5A we see that the series X^^g M*[p(|x|)^ + 2, 1] must con- 
verge, since the sum is — /? plus a sum over probabilities of mutually exclusive events. Furthermore, 
we have 

Pr[3 1 : = • • • = Rx^t-i = 0, Rx,t = !]>/? 
if and only if J2t>o + 2, 1] > 0. Theorem |]l now follows from Theorem p. 6 
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